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Periodic orbit action correlations are studied for the piecewise linear, area- 
preserving Baker map. Semiclassical periodic orbit formulae together with universal 
spectral statistics in the corresponding quantum Baker map suggest the existence 
of universal periodic orbit correlations. The calculation of periodic orbit sums for 
the Baker map can be performed with the help of a Perron-Frobenius type operator. 
This makes it possible to study periodic orbit correlations for orbits with period up 
to 500 iterations of the map. Periodic orbit correlations are found to agree quan- 
titatively with the predictions from random matrix theory up to a critical length 
determined by the semiclassical error. Exponentially increasing terms dominate 
the correlations for longer orbits which are due to the violation of unitarity in the 
semiclassical approximation. 



I. INTRODUCTION 

Numerical evidence suggests that eigenvalue spectra of individual low dimensional quantum 
systems show correlations which are solely determined by the underlying classical dynamics, sym- 
metries and Planck's constant h (see e.g. Bohigas et al 1984, Berry 1987, Bohigas 1991). Especially 
quantum systems exhibiting fully chaotic or completely regular dynamics in the classical limit show 
spectral fluctuations coinciding with random matrix theory (RMT). Semiclassical trace formulae 
provide a link between the set of eigenenergies and the set of all periodic orbits in the classical 
system and have proven to be important in understanding universality in spectral statistics and 
deviations thereof (Hannay and Ozorio de Almeida 1984, Berry 1985, Berry 1988, Aurich and 
Steiner 1995, Bogomolny and Keating 1996a). A complete description of spectral properties in 
terms of the underlying classical dynamics of the system is, however, still missing. 

Universal spectral statistics is intimately connected to correlations in periodic orbit length or 
actions distributions. The validity of the random matrix conjecture implies generalised universal 
periodic action correlation for systems for which the trace formula is exact (Argaman et al 1993). 
Keating (1993) applied this idea to the Riemann zeta function whose non-trivial zeros are con- 
jectured to follow random matrix statistics. Here, the Hardy - Littlewood conjecture provides an 
explicit expression for the pair - correlations of prime numbers, (the analogue of periodic orbits 
in the Riemann case). Expressions for the two-point correlation function of Riemann zeros can 
be obtained making use of the conjectured prime number correlations which coincide asymptot- 
ically with RMT results for an ensemble of Gaussian random matrices invariant under unitary 
transformation (GUE). A generalisation of these arguments to n-point correlation functions for 
Riemann-zeros has been given by Bogomolny and Keating (1995, 1996). 

Far less is known about action correlations of periodic orbits for generic chaotic systems and even 
numerical evidence for the existence of these universal periodic orbit correlations is weak. Argaman 
et al (1993) present numerical results which agree qualitatively with the RMT-prediction presented 
in the same paper; Cohen et al (1998) find correlations in the length spectrum of periodic orbits 
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in the Stadium billiard which indicate RMT-like correlations but a quantitative analysis can not 
compete with the accuracy reached for spectral eigenvalue statistics. Other groups report complete 
failure in their search for any kind of action correlations (Aurich 1998) or find action correlations 
which do not coincide with the RMT - prediction by Argaman et al (1993), see Sano (1999); in 
addition there are yet no arguments other than from the duality between quantum eigenvalues and 
classical periodic orbits which indicate why universal classical correlations should exist. Cohen et 
al (1998) present a variety of ideas based on identifying relationships within various subgroups of 
orbits in the Sinai billiard. The results are, however, still of speculative nature. 

In this paper I will study quantum and scmiclassical spectra and discuss necessary conditions 
for the existence of universal periodic orbit action correlations. I will focus in particular on the 
classical and quantum Baker map. It will become evident that the periodic orbit correlations 
deviate from the universal behavior for long orbits due to the violation of quantum unitarity in 
the scmiclassical approximations. Correlations following the RMT-prediction do, however, exist 
for short periodic orbits which can not be explained by classical sum rules alone. Properties of 
periodic orbits can be studied in detail for the Baker map by writing periodic orbit sums in terms of 
a suitable classical Perron-Frobenius operator. Exponentially increasing terms in the large action 
limit can be regularised by imposing unitarity onto a scmiclassical quantization as proposed by 
Bogomolny and Keating (1996a). 



II. PERIODIC ORBIT CORRELATION FUNCTIONS FOR QUANTUM MAPS 

In the following I will limit the discussion to quantum maps. A quantum map acts on a finite 
dimensional Hilbert space and the quantum dynamics is governed by the equation 

Ipn+l = (1) 

with U being a unitary matrix of dimension N . We assume that the map has a well defined classical 
limit for N oo described by a discrete dynamical map. ipn is the discretised N - dimensional 
wave vector and N is equivalent to the inverse of Planck's constant h, i.e. N 1/h. For possible 
quantization procedures of classical maps see e.g. Hannay and Berry (1980) for the cat map, Balazs 
and Voros (1989) and Saraceno and Voros (1994) for the Baker map or Bogomolny (1992), Doron 
and Smilansky (1991), and Prosen (1994, 1995) for scmiclassical and quantum Poincare maps. The 
spectral density of eigenphases 0i of U can be written in the form 



N ^ ^ oo 

d{e, N)^y] 52.{e -e,)^ — + -Re V TrU"e-™^ (2) 

i—1 n—1 

with d — N/2ti being the mean density of eigenphases in the interval [0, 2n\ and 52-r: denotes the 
periodically continued delta function with period 27r. The spectral measure of interest here is the 
two-point correlation function defined as 

i?2(a;, ^) = ^ < d{e)d{e + x/d) > . (3) 
d 

The average is taken over the unit circle, i.e. < . >= J^^ . dO which leads to 

R2{x, - - ^) = 1 + ]^ E |TrU"(7V)|2 cos(27r^a;) ; (4) 

— l n—1 

the spectrum {xi — dOi] i — 1,2... N} is unfolded to mean level density one. Note that i?2, as 
defined in (|^), is still a distribution and further averaging has to be applied. Fourier transformation 
of the two-point correlation function yields the so-called form factor which can be written as 
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the variable r takes the discrete values r — n/N, rt, e IN which reflects the periodicity of the 
two-point correlation function R2. Note that the form factor does not converge in the limit N ^ 00 
but fluctuates around a mean value. The limit distribution of K{t, N) ioi N ^ 00 is, however, 
expected to exist and is conjectured to be Gaussian (Prange 1997). I will refer to the form factor 
as the mean value of the distribution near a given t value, i.e. further averaging over small 
T-intervals has to be performed. 

The traces of U" can in semiclassical approximation be written as sum over periodic orbits of 
topological length n of the underlying classical map (Gutzwiller 1990, Bogomolny 1992), i.e. 

(n) 

TrU"(iV) w TrU"(7V) := ^ Ape^"^^" . (6) 

p 

The complex prefactors Ap contain information about the stability of the periodic orbit and Sp 
denotes the action or generating function of the map corresponding to the periodic orbit p. The 
relation between quantum traces and semiclassical periodic orbit sums is in general obtained by 
stationary phase approximation and is exact only in the limit t^^ = N/n —^ 00 (Andersson and 
Melrose 1997). Exceptions thereof are e.g. the cat map (Hannay and Berry 1980, Keating 1991) 
or quantum graphs (Kottos and Smilansky 1997) for which Eq. (Q) is exact for all r. 

The relation (||) makes it possible to connect purely classical action correlation functions with 
the statistics of quantum eigenvalues. A weighted action correlation function can be obtained by 
considering 

00 (n) 

P{s,n)^ |TrU"(iV)|V-^^^ = ^Ap^;,,52.(5p-^p. -s) (7) 

N——00 p^p' 

which is a Fourier sum in N depending only on classical quantities such as the topological length 
of the orbit, n, and the actions Sp and weights Ap. Eqs. (^) and (0) together with indicate 
that energy level statistics and periodic orbit action correlations are connected provided that both 
the quantum traces TrU" and the semiclassical 'traces' TrU" follow the same probability distribu- 
tions, i.e. random matrix statistics. Little is known about the statistical properties of semiclassical 
traces and eigenvalues and I will study the existence of universal RMT-behavior for semiclassical 



expressions in more detail in section IV 



The RMT-conjecture implies for quantum traces TrU"(A^) of systems whose classical limit is 
chaotic and non-time reversal symmetric 

lin, 1 < |TrU"(iV)p >= ( " ^ ^ ; (8) 

jv-,00 A'^ ' [I for n > A ' ^ ' 

the brackets < . > indicate the average over an n-interval small compared to the dimension N. 
Assuming that the above relation holds also for the semiclassical traces < |TrU"(A^)p > we can 
insert (||) in (|^) which leads to the asymptotic result (Argaman et al 1993) 

N ^ / sinn7rs\^ ^, ^ , , , ^ 

P{s,n) = P{n)-{——\ +Ji5{s) for |s| < 1. (9) 

Here, P{n) — X)p"p' -^p-^p' the mean part of the weighted periodic orbit action pair density. 
Note that the number of periodic orbits increase exponentially with n for chaotic maps which in 
turn implies an exponential increase for the density of periodic orbit actions (modulus 1) and in 
general also for the weighted density P{n). 

After rescaling the s - variable according to s ^ a/n, Eq. (|^) can be written as 



- (-) - ^ - ( ^ ) + Sia) for |.| « n, (10) 



2 

which puts the non-trivial, cr-dependent correlations into n - independent form 



3 



The corresponding equation for systems with time-reversal symmetry is (Argaman et al 1993) 



2 

H [cos 27rcr (si(27rf7) cos 2TTa — Ci(27rt7) sin 2TTa) 

TTCT 

+ Ci(47rcr) sin47rCT - si(47rcr) cos47rf7] + 2S{(t) for \a\ < n, (11) 

with si (a;), Ci(x) being the Sine - and Cosine - integrals. 

For integrable systems conjectured to have Poissonian spectral statistics one expects 

(") 

— < |TrU"(iV)P >=< ^ApA;,e2"^(^^-^^') >= 1 for \N\ ^ 
p,p' 

and thus 

oo ^ 

P^°'^(s,n)= ^|TrU"(7V)|2e2'^*^^ =P(n)-l + ,5(s), (12) 

N=-oo 

i.e. periodic orbit actions are uncorrelated. 

The (5-functions in (^ - (|l^) can be identified with the diagonal terms p — p' in the sum (|^) 
after making use of the classical sum rule (Hannay and Ozorio de Almeida 1984) 

(n) 

^lApp^l for large n. (13) 

p 

Before investigating the existence of correlations of the form (|l^) or (111]) in more detail, I 
would like to make a few remarks: the weighted periodic orbit correlation function contains a 
background term P{n) which increases in general exponentially with n. This is in contrast to the 
linear behaviour of the mean quantum level density d{N) in (||) and (^. Note also, that the complex 
weights Ap carry phases which may result in cancellations in P and may occasionally lead to P = 
0. The correlation functions (|l^) and (|l|) are in addition not rescaled with respect to the mean 
action density as e.g the two-point spectral correlation function i?2 in (^). The differences in action 
between adjacent orbits of the same length n is exponentially small on the s or ct scale and the 
term — (sinTrcr/Trcr)^ in (p^, ( pi] ) indicates long range correlations on the scale of the mean periodic 
action separation. The correlations ( p^ ) and ( pi] ) imposed by random matrix theory are thus a 
small modulation of order 0(1) on top of an exponentially large background when considering all 
periodic orbit pairs (Dahlqvist 1995). The two-point correlation function for periodic orbit actions 
approaches the Poissonian limit for n — > oo when rescaling P{a, n) with respect to the mean action 
density in accordance to e.g. the definition of the spectral correlation function (^); periodic orbits 
are thus uncorrelated on scales of the mean action separation in agreement with the numerical 
results by Harayama and Shudo (1992). 

The search for an C(l) - effect on top of an exponentially large background is one of the 
main obstacles when investigating RMT-induced classical action correlations numerically. The 
second major problem in studying periodic orbit properties especially in the large n - limit is the 
exponential increase in the number of orbits with increasing n. In the following I will focus on a 
specific example, the classical and quantum Baker map, in which the problems mentioned earlier 
can be overcome by constructing a quasiclassical Perron-Frobenius type operator (Dittes et al 
1994). 

III. THE BAKER MAP 

The Baker map has become a standard example when studying classical and quantum chaos (see 
e.g. Balazs and Voros 1989, Saraceno and Voros 1994, O'Connor et al 1992, Hannay et al 1994). 
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The classical dynamics is given by a two-dimensional, piecewise linear, area-preserving map on the 
unit square defined as 



q' = 2q-e 

p' = i(p + e) with e=[2q] (14) 

and {q,p), {q' ,p') being the initial and final points. The notation [x] stands for the integer part of x. 
The Baker map is hyperbolic with a well defined Markov - partition and a complete binary symbolic 
dynamics given in terms of the e = {0, 1}. The map is invariant under the anti-unitary symmetry 
T : {q,p) —^ (j>,q) (being equivalent to time reversal symmetry) and a parity transformation 
i? : ((7,p) — > (1 — g, 1 — p). A proper desymmetrisation of quantum spectra and semiclassical 
periodic orbit sums with respect to parity will be crucial when studying spectral statistics and 



periodic orbit correlations in section IV 



Each periodic orbit of the map can be associated with a finite symbol string e = (ei, £2, . . . , £«) 
and the phase space coordinates of the orbit are given by the expression 



9^ = F^E^^2— 1; p,._L_g,,2V (15) 

1=0 1=0 

A suitable generating function of the map is 

W,{p',q) = 2p'q-ep' -eq with e = [2g] , (16) 

which in turn defines the action; the action of a periodic orbit e of the map can up to an additive 
integer be written as (Dittes et al 1994) 



l)e, = V(g,-l)[2<z,] (17) 



i=l i=l 



with qi being the g-coordinates of the periodic orbit. The sum in ( pT[ ) is invariant under time 
reversal symmetry and parity-transformation R{qi) = \ — qi\ the later can be shown using the 
relation Ya^-^ qi = Ya^i Ci- 

A quantized version of the Baker map is obtained by making use of a discretised version of the 
generating function (|l^) and the classical structure of the map in the mixed representation. A 
suitable formulation preserving all the classical symmetries is provided by the choice (Balazs and 
Voros 1989, Saraceno and Voros 1994) 



V/2 



U(iV) = F-x( ^J/^ ) (18) 



(F^)ij = 1 e-2"(^-^)a-^)/^ i, j = 1, . . . . (19) 

V A 



with 



The Fourier-matrix F is nothing but the transformation from position to momentum representation 
in the finite dimensional Hilbert space. The dimension A^ of the map U(A^) has to be even in this 
construction and A^ is equivalent to the inverse of Planck's constant. The unitary map U(A^) 
commutes with the parity operator (R7v)i,j = Si^N-j+i due to the particular choice of half-integer 
phases (Saraceno and Voros 1994). Symmetry reduction is obtained by considering the matrices 
U±(A^) defined as 

U± (A^) = i (U(7V) ± U( A)R^) (20) 

acting on the symmetric/anti-symmetric wave- vectors only. This leads effectively to a reduction 
in dimension by a factor of 2. The traces of U!J- can in the large A^-limit be written as (Saraceno 
and Voros 1994, Toscano et al 1997) 
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1 (") / 9"/2 9n/2 



e 



(21) 

The first two terms in tlie sum over the possible finite symbol strings of length n are the usual 
semiclassical Gutzwiller-like periodic orbit contributions. They arise as stationary phase points in 
a continuum approximation of TrU", Tr(U"R), respectively, i.e. sums over matrix elements are 
replaced by integrals (Tabor 1983). The action Se of a periodic orbit of length n is given by Eq. (pT|). 
S'e corresponds to half the action of an orbit of length 2n with symbol code (ei, . . . e„, ei, . . . £„) 
and ei = 1 — e^. These are the orbits being invariant under the classical parity transformation 
R and are hence the stationary phase contributions of Tr(U"R). The anomalous aglogiV + be 
corrections are due to the discretisation of the phase-space and arise in 'sub leading' integrals when 
performing Poisson-summation on TrU" i.e. replacing sums by sums over integrals (Saraceno and 
Voros 1994, Toscano et al 1997)). Next leading terms are diffraction corrections arising from the 



discontinuous nature of the classical and quantum map (14) and (|18|). Note the extra phase tt in 
the Tr(U"R) - contributions which originates from the discretised stationary phase approximation. 

In the following I will mainly concentrate on the Gutzwiller - like periodic orbit contributions, 
i.e. I will consider the approximation 

on/2 (") 

TrU" (iV) « TrU" {N) = — ^ e^'^'^^e (22) 



TrU!^(iV) « TrU!^(iV) = i 



e 

on/2 (") 9n/2 (") 

2n_lZ^ "'"2" + l^ 



I will henceforth refer to the periodic orbit sums (p^ ) as the semiclassical approximation of the 
quantum Baker map. The influence of the remaining contributions in (|2l|), including the leading 
order (log N) - corrections, will be discussed at the end of this section. 

The main problem in calculating periodic orbit sums like ( p2| ) is the exponential increase of 
periodic orbit contributions. Even for the Baker map where periodic orbit actions are given by the 
simple analytic formula ( [l7| ) a direct summation of TrU" becomes a computational challenge for 
n > 30. For the Baker map it is, however, possible, to construct a quasiclassical operator U whose 
traces coincide with the semiclassical periodic orbit sum ( p^ ) (Dittes et al 1994). The spectrum of 
this operator can be computed explicitly and simplifies the task of calculating traces, i.e. periodic 
orbit sums, considerably. Other quasiclassical operators for the Baker map have been proposed 
by Kaplan and Heller (1996), Sano (1999) and Hannay (1999). Note, that an identification of 
semiclassical periodic orbit sums with classical or quasiclassical operator is not possible for generic 
systems and the notation TrU" is then a mere substitution for the periodic orbit sum itself. 

The Baker map is special for the following reasons: firstly, the periodic orbit amplitudes Ap, 
(cf. Eq. (||)), do not depend on the specific orbit, but only on the topological length n, which is 
a consequence of the piecewise linearity of the classical map. Secondly, the dynamics of the q - 
coordinate as well as the periodic orbit - action ( pT| ) are independent of the momentum p which 
allow for a classical separation of momentum and configuration space variables. As a consequence, 
one can define a one dimensional Perron-Frobenius type integral kernel acting on the g-coordinate 
only (Dittes et al 1994). In a form preserving the parity-symmetry R{q) = 1 — q it can be written 
as 

U{q, q'; N) = V2 S {q - {2q' - [2q'])) e^-^^^C?') (38) 

with 

S{q')=q'[2q']-^i[2q']+q'). (24) 
One easily deduces U{q,q') — U{Rq,Rq'). The traces of IJ{q,q') coincide with the periodic orbit 



sum (22) if the operator is defined on the space of analytic functions on the unit interval with peri- 



odic boundary conditions (Dittes et al 1994, Rugh 1992). The operator is thus infinite dimensional 
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(having infinitely many eigenvalues) in contrast to the finite dimensional quantum matrix U(iV) 
and A'^ is a mere parameter in (^3|). A symmetry- reduced version of this quasiclassical operator is 
obtained by considering 



U±{q. q'; N) = ^ (Uiq, q'; N) T U{q, Rq'; iV)) 



(25) 



Note, that the quasiclassical operator acting on antisymmetric functions (denoted U+ here) cor- 
responds to the symmetric quantum subspace U+ and vice versa, reflecting the 'quantum' origin 
of the extra phase for the periodic orbit contributions to TrU"R, see Eq. (pi|). The operator (|2^) 
or ( p5| ) can now be represented as an infinite dimensional matrix after choosing a suitable basis 
set and the calculation of traces can be reduces to matrix calculus (provided the operator is trace 
class). An obvious choice for the basis functions is the Fourier-basis, i.e. 



Uk.n= / dq d(7';7(g,q';iV)e2-("«'-'=«) 

Jo JQ 



(-1)'"V2- 



-TriN/2 
2TTi 



N/2 + m-2k 



N/2 -m + 2k 



Similarly one obtains 



{UR)k,rn= dq dg't/(g,i?g';iV)e2-('"«'-'=«) 
Jo Jo 



-l)"\/2. 



,-niN/2 

27ri 



„7ri(JV/2+m+2*;) _ ^ ^TTi{N /2-m-2k) _ ^ 



N/2- 



2k 



N/2 -m-2k 



(26) 



(27) 



with fc,m integers. The various terms (e^*(^/2±™±2fe) _ iy(^]\[^2 ± m ± 2k) are equal to it: for 
N/2±m±2k = 0. 



a) 



ImA 




b) 



ImA 




ReA ReA 

FIG. 1. Semiclassical (o) and quantum {+) eigenvalues of the Baker map for A*' — 
subspaces are shown separately: a) positive parity, b) negative parity. 



80; the parity 



A typical eigenvalue spectrum {A.;}, i = 0, 1, 2 . . . of the quasiclassical operator U± is shown in 
Fig. |l| together with the quantum eigenvalues. In the following, the eigenvalues will be ordered 
with decreasing modulus, i.e. \Ai\ > \Aj\ for i < j. The operator U approximates unitarity of the 
corresponding quantum map U quite well; approximately eigenvalues in each parity subspace 
lie near the unit circle and these eigenvalues agree well with quantum results. The modulus of 
semiclassical eigenvalues becomes exponentially small for i > N/2 in each subspace and infinitely 
many eigenvalues 'disappear' by spiraling into the origin. This behavior reflects the trace-class 
property of the operator and makes it possible to truncate the size of the matrix U± in actual 
calculations. (Choosing dim(U±) = 3 x ^ ensures convergence of the first N/2 eigenvalues in each 
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parity subspace to 4 significant digits). 



0.14 
0.12 R 

0.1 d ~ n' 

0.08 * 
0.06 i: 
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+ pos. parity 
□ neg. parity 




50 



150 
N/2 



250 



300 



FIG. 2. The semiclassical eigenvalue with the largest modulus plotted here as 7*(A'^) = log|Aj(A'^)| 
versus N/2 for both subspaces; the deviation from the unit circle decreases in both cases like l/y/W 
(dashed lines), but one finds typically 7+ (A'') < j^lN). 

A measure for the error of the semiclassical approximation compared to quantum calculations 
is the deviation of the semiclassical eigenvalue with largest modulus, Ao(A^), from the unit circle. 
The distribution of ^^{N) := log |AJ(A/)| for integer values of N and both subspaces is shown in 
Fig. ||. The deviation of semiclassical eigenvalues from the unit - circle decreases with increasing N 
following a 1/\/N behavior. The - scaling suggest that corrections to individual semiclassical 

eigenvalues are dominated by diffraction contributions. The log(iV) corrections in the traces, 
cf. (21), become dominant when summing over the periodic orbit contributions or equivalently 
over the semiclassical eigenvalues and are thus a collective effect of the various contributions. 

The deviations from the unit circle for semiclassical eigenvalues with modulus greater than one 
is systematically larger for eigenvalues corresponding to the C/_ - operator compared to those from 
the [/-!_ - operator. This can qualitatively be understood for small N when considering the limit 
= 0; the operator U{q,q'\N = 0) = \/2S{q — {2q' — [2q'])) is the classical Perron-Frobenius 
operator for the Sawtooth-map up to a factor of \/2. The spectrum U{q, q'; 0) has therefor a single 
non-zero eigenvalue V2 in the symmetric subspace which corresponds to U- in our notation all 
other eigenvalues are zero. The eigenvalues are smooth functions of the continuous variable N 
which leads to |Aq | > \ for A^ smaller one. The fact that this behavior persists for large A^ 
values , cf. Fig. ^ can certainly not be explained by the argument above but is worth noting as a 
purely numeric result. 



IV. SEMICLASSICAL SPECTRAL STATISTICS AND ACTION CORRELATION 
FUNCTIONS: NUMERICAL RESULTS 

The existence of a quasiclassical operator for the Baker map makes it possible to write traces as 
convergent sums over the quasiclassical eigenvalue spectrum, i.e. 

00 

TTUi{N)^Y.i^tWr ■ (28) 

1=0 

The violation of unitarity and the existence of semiclassical eigenvalues with modulus greater than 
one leads to the asymptotic result 

\TtU!^{N)\ ^ exp{n-/^{N)) for n -> 00 
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with 7^(A^) — log|Aj(A^)|, i.e. the traces grow exponentiaUy in the hmit n ^ oo (Keating 1994, 
Aurich and Sieber 1994). This kind of behavior is connected to the notorious convergence problems 
for periodic orbit sums of the form 

oo (n) 
n— 1 p 

which are absolutely convergent only for complex N with Ini{N) > ht/2 and ht is the topological 
entropy (Eckhard and Aurell 1989). For the Baker map the absolute sum over periodic orbit 
contributions corresponds to Trf/"(0); this sum diverges like 2"/^ for large n, i.e. ht = 

log 2 here. A summation over the semiclassical eigenvalues (^8|) gives a detailed account of the 
regime of conditional convergence for periodic orbit sums as a function of the real part of N. 

The exponential increase in the traces ( |28| ) leads to exponentially increasing terms in a semi- 
classical approximation of the form factor 



Ksc{r,N) = ^iTrJ/^p = ^ E {^t^Y^ « exp(27VT7±(A^)) for r ^ oo . (29) 



A stationary distribution of Kgc for fixed r and iV ^ oo is obtained if < j{N) > falls off faster 
than A''"^. Such a behavior is expected for smooth hyperbolic maps where the semiclassical error 
is dominated by 1/A'^ - corrections due to higher order terms in a stationary phase expansion. This 
is, however, not the case for the Baker map where diffraction effects of the order N~^^^ dominate, 
see Fig. || 

Using the semiclassical traces Tr[/" to evaluate the two-point correlation function (^ will, on 
the other hand, always be affected by the exponentially increasing terms in the form factor; a semi- 
classical calculation of R2{x, N) diverges thus independently of N and the behavior of < 7(A^) >. 

The weighted classical correlation function (|7|) is a finite sum and thus well defined for fixed 
topological length n. From (|25) one obtains in the large n - limit 



P{n) — 2 " in the '+' - subspace 
P{n) — 2" in the '— ' - subspace 

and we will discard this 'trivial' mean part from now on. The large n - asymptotics of the non- 
trivial correlations in (^ is dominated by the eigenvalue being the largest among all Ao{N); after 
defining 

7o± :=log|A±(7Vo)| := sup (log |A±(A^)|) 

where A^o is the A^ - value corresponding to the largest eigenvalue Aq, one obtains 

P(s, n) ^ e^"^"* cos(27riVos) ^ oo for n ^ oo . (30) 

It follows from the considerations above that a full description of quantum spectral statistics 
in terms of classical actions and amplitudes needs to incorporate unitarity into a semiclassical 
approximation (Keating 1994). I will come back to this point at the end of the section. Apart 
from that there seems to be little information in the asymptotic behavior of Ksc(t, N) or P{s, n). 
The rates |7*(Af)| > correspond to 'semiclassical escape rates' and contain little information 
about the classical dynamics or the quantum map. The distribution of these rates is instead of 
statistical nature due to accumulation of errors in the semiclassical approximation. 

The more interesting question in this context is whether or not semiclassical periodic orbit 
formulae for classical chaotic systems are able to reproduce RMT-statistics quantitatively in the 
non-asymptotic regime. The form factor obtained from the quasiclassical operator is displayed 
in Fig. H together with the result obtained from the quantum spectrum and the GOE-prediction 
(applicable for systems with time reversal symmetry). The semiclassical result is indeed capable of 
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reproducing the RMT-result in both subspaces for smaU r values; especiahy the deviation of the 
RMT-behavior from the classical result K(t) — 2t obtained from the sum rule (|l^) in the limit 
T ^ is here well reproduced by periodic orbit formulae. This clearly indicates that there are 
non-trivial correlations between periodic orbit actions for systems with time-reversal symmetry. 




b) 



K 




1 1.5 2 2.5 0.5 1 1.5 2 2.5 

X = n/N T = n/N 

FIG. 3. The form factor K{t), i.e. the square modulus of the quantum traces TrU" (full line) and the 

semiclassical traces Trt/" (dashed line) for positive (a) and negative (b) parity is displayed as function of 
T = n/N for N/2 = 259 together with the RMT - result for the Gaussian Orthogonal Ensemble (smooth 
curve). The semiclassical break-times, Eq. (^), are ~ 1.002, t/^ ~ .301, respectively. 

Representing periodic orbit formulae by sums over quasiclassical eigenvalues allows one to study 
the large N and r limit of a semiclassical approximation of the form factor. Exponentially increas- 
ing components start to dominate Ksc{t, N) for t larger than a semiclassical break-time of the 
order 



1 



7±(iV)iV 



The break-time is typically a factor 5 larger in the 
Fig. |. 



subspace than in the '- 



(31) 

subspace, see 



The action correlation function (j^), on the other hand, can be obtained either by sampling 
action-differences directly (which is limited to n < 15-20 due to the exponential increase in the 
number of orbits) or with the help of the traces Tr[/" [N) . The latter method demands calculation 
of the spectra of fj (N) for integer A^- values and summation of the Fourier-sum in (|^) directly. In 
practice, the sum is truncated at a finite iV-value, iVmax, which corresponds to a smoothing of the 
original correlation function P{s,n) on scales l/A^max; after rescaling according to Eq. ( [lO| ) and 
subtracting P{n) one obtains 



n ^ina x AT 

Pscai{a,n,Nn,.^) = — y |TrU"(Ar)|2cos(27r-a) 



with 



sin(27ra2:) 



and a — 



da'Pscai{cT',n)gaicr- a') (32) 



(33) 



The smoothed periodic orbit correlation function Pscaii<^, n, A'max) is shown in Fig. ^for various 
n-values and A'max = n. (Contributions from diagonal terms leading to the S - functions in (|^) 
and ( p^ are subtracted here). Universal periodic orbit correlations in the subspace is observed 
up to n « 70; these are 2"^° distinct periodic orbits, a number inaccessible to pure periodic orbit 
calculations. One can study even longer orbits and finds exponentially increasing terms dominating 
the correlation function for n values above the transition point n « 70. The modulations in the 
periodic orbit pair density function P{cr, n, n) are indeed orders of magnitudes larger than the 
RMT-correlations for the 2^°" periodic orbits of length n = 500, see Fig. ^ 
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Things look similar in the '— ' subspace, the deviation from universality occur, however, for much 
smaller n-values, i.e. at n w 6. The correlations which occur for larger n have the simple form 
P(s, n) — e^o " cos(27rs), see Fig. for n = 65. This oscillatory behavior has already been observed 
by Sano (1999) when studying periodic orbit correlations in the Baker map without separating the 
two symmetry subspaces. Sano's result can now be interpreted in terms of the spectrum of the 
quasiclassical operator ([23|); the exponents 7^ for the two different subspaces, which determine 
the asymptotic behavior of P(s, n) for large n, see Eq. (|30|), are 

7fJ" = 0.022885 with Nq = 14 
7(7 = 0.148508 with iVo = 1 , 

which can also be deduced from Fig. ^. The exponentially growing terms in the periodic orbit 
correlation function have thus a seven times larger leading exponent in the '— ' subspace compared 
to the '+' subspace. The correlation function P{s,n) in the '— ' subspace is dominated by the 
Fourier coefficient |Tr?7"'(l)p for n > 6 and RMT-type correlation can not develop. The same is 
obviously true when considering the full Baker map without symmetry-reduction. 






CT a 
FIG. 4. The periodic orbit correlation function are shown for various n - values in the '+' subspace 

(Fig. a, c) and ' — ' subspace (Fig. b, d); 'universal' periodic orbit correlation which coincide with the GOE 

- prediction can be observed up to n « 70 in the '+' subspace, but only up to n « 8 in the ' — ' subspace. 

Non- universal periodic orbit correlations develop for large n (with exponentially increasing amplitude), see 

Fig. c and d. The cut-off is A'^max = n here, see Eq. (p2). 



Universal correlations coinciding with the RMT-prediction do, however, exist. This is clearly 
demonstrated in the subspace for n < 70. It is the violation of unitarity in a semiclassical 
approximation which leads (in general) to exponentially diverging terms which in turn overwhelm 
the pair-distribution function P{s, n) for large n. The existence of universal periodic orbit corre- 
lations is thus strongly linked to the preservation of unitarity in a semiclassical formulation. The 
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absence of RMT - like periodic orbit correlation, as observed in the ' — ' subspace, on the other 
hand, contains httle information about the phenomena of universal spectral statistics but merely 
reflects the limitations of the approach due to the underlying semiclassical approximations. 

A possibility to enforce unitarity onto a semiclassical description has been proposed by Bo- 
gomolny and Keating (1996a). The starting point is the quantum staircase function Af{9) = 

^f^i QiO — O.j) with Qi being the quantum eigenphases and Q{x) denotes the Heaviside step- 
function. A smoothed, semiclassical version of the staircase function may be written as the trun- 
cated Fourier-series of M{9) 

N 1 ^1 
= — - -/m^ -TrU"(Ar)e-"^ , (34) 

and the traces TrU" are expressed as sums over periodic orbit contributions, see Eq. (H). The 
unitarity condition is implemented by choosing a quantization condition 

AA(^„,7V) =n+i (35) 

with On being the solutions of ( ^ ) ; the On represent a semiclassical approximation of the quantum 
eigenphases On- A density of states is defined according to 

N ^ / 1\ 

D{0, N)^J2 - = diO, N) J2 52. {M{0, N)~n^-\ (36) 
1=1 1=1 ^ ^ 



with d{0,N) = -^^[{9, N). The density ( p6[ ) is again written in terms of periodic orbit contri- 
butions only but 'preserves' unitarity by construction. Bogomolny and Keating (1996a) started 
from the expression ( p^ ) to derive periodic orbit corrections to the quantum two-point correlation 
function i?2 beyond the results obtained from the classical sum-rule (p^. 

The Perron-Frobenius representation of the semiclassical traces for the Baker map allows one 
to test the quantization condition (^5|) for large up to iV w 500 by calculating the semiclassical 
traces in (53) directly; the 'regularised' spectrum {Oi} thus obtained is used to construct new 

— ~ n 

semiclassical traces TrU = '^^exp(in9i) which in turn are used to calculate the periodic orbit 
correlation function (^. Note, however, that the simple relation between periodic orbits of length 
n and the n-th semiclassical trace gets lost here and the new periodic orbit correlation function 
P(s, n) contains also correlation between orbits and pseudo - orbits (being composites of shorter 
orbits). 

The correlation function obtained from the 'regularised' semiclassical data is shown in Fig. |^ 
together with correlation functions where regularisation has not been applied, cf. also Fig. ^. The 
functions P{(J, n, A^max) are shown here with enhanced resolution compared to Fig. ^, i.e. with a 
cut-off A^max — ^n. The periodic orbit correlation function for n — 70 follows the random matrix 
prediction also with enhanced resolution, but again deviates from the universal behavior for n = 
200. The exponentially growing terms are eliminated in the corresponding regularised correlation 
function using the spectrum obtained from (|35| ) for n = 200. The regularisation procedure unveils 
the underlying universal correlations in the classical actions, see Fig. ^. 

The regularisation process ( ^5|) has the disadvantage that the Fourier-components of the new 
density of states D{0, N) can no longer be written as closed expressions of a finite number of periodic 
orbits; they become complicated infinite periodic orbit sums after writing the 827^ - function in 
its Fourier-components (Bogomolny and Keating 1996a). The resulting periodic orbit correlation 
functions, as e.g. shown in Fig. |^ for n = 200, consist of periodic orbit contributions of orbits and 
pseudo - orbits of topological length up to and including n = 200 which does not simplify the task 
of understanding the origin of universal periodic orbit correlations. 
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FIG. 5. The periodic orbit correlation function shown for various n - values in the '+' subspace with 
enhanced resolution compared to Fig. ^; the regularisation procedure eliminates the exponentially 
growing terms and universal correlations are uncovered now also for n = 200. 

V. CONCLUSIONS AND OUTLOOK 

Correlations between actions of periodic orbits of the Baker map up to orbits of topological 
length n = 500 (corresponding to 2^"'' different periodic orbits) have been studied with the help 
of a quasiclassical Perron-Frobenius operator. The spectral form factor can be calculated by 
purely semiclassical expressions which coincide for small r with random matrix results beyond 
the validity of the classical sum-rule, but diverge for t oo. Action correlations of periodic 
orbits have been investigated which show universal non-trivial correlations linked to random matrix 
theory for short periodic orbits, but depart from the universal behavior for long orbits due to the 
violation of unitarity in a semiclassical approximation. The transition point from universal to non- 
universal statistics is distinctively different for the two parity subspaces in the Baker map. It is 
linked to the magnitude of the semiclassical error made when approximating the quantum density 
of states by semiclassical periodic orbit formulae and is controlled by the largest deviation of a 
semiclassical eigenvalue from the unit circle. This behaviour is probably generic for quantum maps. 
Statistical properties of semiclassical expressions are universal in the non-asymptotic regime, the 
transition point at which universality breaks down is, however, system dependent and controlled 
by the semiclassical error. Imposing unitarity on a semiclassical approximation makes it possible 
to discard exponentially growing terms in the periodic orbit correlation function which in turn 
uncovers universal correlation even above the transition point, as demonstrated here for the Baker 
map. 

The study presented here establishes for the first time unambiguously the existence of universal 
periodic orbit correlations in a classical chaotic system whose quantum counterpart shows RMT- 
eigenvalue statistics; the limitations due to the semiclassical approximations are discussed in detail. 
This sets a proper framework for studying the origin of these classical correlations. A better under- 
standing of the interplay between classical periodic orbit action correlations and unitarity might 
shed light on the existence of universality in quantum spectra in general. 
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